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Abstract. We consider searches for dark matter annihilation to monoenergetic neutrinos 
in the core of the Sun. We find that liquid scintillation neutrino detectors have enhanced 
sensitivity to this class of dark matter models, due to the energy and angular resolution 
possible for electron neutrinos and antineutrinos that scatter via charged-current interactions. 
In particular we find that KamLAND, utilizing existing data, could provide better sensitivity 
to such models than any current direct detection experiment for mx 15 GeV. KamLAND’s 
sensitivity is signal-limited, and future liquid scintillation or liquid argon detectors with 
similar energy and angular resolution, but with larger exposure, will provide significantly 
better sensitivity. These detectors may be particularly powerful probes of dark matter with 
mass 0(10) GeV. 
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1 Introduction 

If the dark matter-nucleon scattering matrix element is spin-dependent, then one of the lead¬ 
ing detection strategies is to search for the neutrinos that arise from the annihilation of dark 
matter that has collected in the Sun [1-3]. Some types of neutrino detectors, such as liquid 
scintillation or liquid argon detectors, can provide very good energy resolution for an elec¬ 
tron (anti)neutrino. In this work, we describe the enhanced sensitivity of these detectors to 
dark matter candidates that annihilate to monoenergetic neutrinos, and in particular deter¬ 
mine the current limits that can be placed by the Kamioka Liquid Scintillator Antineutrino 
Detector (KamLAND) with data that has already been taken. 

WIMP annihilation directly to neutrino-antineutrino final states, XX —> vv, is typically 
considered to be disfavored. Indeed, if the dark matter candidate is a Majorana fermion, 
as is found in the MSSM for neutralino dark matter, direct annihilation to vv is strongly 
suppressed: The initial state consists of two identical fermions, so its wavefunction must be 
totally antisymmetric, implying that s-wave annihilation can only occur from a J = 0 initial 
state. For a J = 0 vv final state, the neutrino and antineutrino must have the same helicity, 
so one would expect this amplitude to be suppressed by the neutrino mass, assuming minimal 
flavor violation. This result relies on the assumption that dark matter is its own anti-particle 
and also on the assumption of minimal flavor violation; if either assumption fails, the process 
XX —> uu may, in fact, have a large branching fraction. Explicit models for which one would 
expect non-negligible annihilation to monoenergetic neutrino final states have been studied 
in [4-9]. 

Searches for neutrinos from dark matter annihilation, in the Sun or elsewhere, are 
sensitive to the expected neutrino energy spectrum. If dark matter annihilates primarily to 
heavy Standard Model particles, such as VF-bosons, 6-quarks, or r-leptons, then the decay 
of these products produces a continuum neutrino spectrum over a range of energies, with 
a high energy cut-off near the mass of the dark matter particle. In this case, experimental 
searches must consider possible excesses of events over the relevant, potentially very broad, 
range of energies. If, on the other hand, dark matter annihilates directly to iw final states, 
the signal would manifest as monoenergetic neutrinos and antineutrinos. 
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The energy of an electron (anti)neutrino can be measured to high precision by liquid 
scintillation detectors [10, 11] and liquid argon detectors [12]. A high energy electron neutrino 
produces an e± through a deep inelastic charged-current interaction, and the energy of the 
electron shower as well as that of the scattered parton can be well-measured at those types 
of detectors. The atmospheric neutrino background to a search for monoenergetic lines is 
greatly reduced as the energy resolution of the detector is improved. Furthermore, if there are 
a sufficiently large number of background events, it may also be possible to characterize and 
subtract the atmospheric background using sideband techniques. Thus liquid scintillation 
or liquid argon detectors may provide impressive sensitivity to dark matter candidates that 
annihilate to monoenergetic neutrinos. 

Previous studies of monoenergetic neutrino signals from dark matter annihilation in 
the Sun have focused on heavier WIMPs accessible to the IceCube neutrino detector. By 
contrast, the detectors considered here are limited in their sensitivity to low WIMP masses 
only by the physics of WIMP interactions within the Sun itself (evaporation) [13], and not 
by a loss of sensitivity to the relevant range of neutrino energies. The detectors considered 
here therefore represent an excellent complement to neutrino searches for heavier WIMPs 
with masses > 100 GeV. 

In section II, we describe the general analysis underlying a monoenergetic neutrino 
search at a liquid scintillation-type detector. In section III, we implement this analysis for 
the specific case of KamLAND. We conclude with a discussion of our results in section IV. 


2 Analysis 

We focus on a search for monoenergetic electron (anti)neutrinos arising from dark matter 
annihilation, XX —> I'li'i, and assume the annihilation cross section is flavor-independent. A 
deep inelastic charged-current interaction (for example, v e d —>• e~u) within the target will 
produce electromagnetic and hadronic showers, the energies of which can be measured quite 
well in liquid scintillation or liquid argon detectors. Moreover, the direction of the electro¬ 
magnetic shower can also be reconstructed 1 , providing a good measurement of the outgoing 
electron or positron’s direction. The energy resolution and angular resolution possible at 
liquid scintillation-type neutrino detectors make these instruments extremely powerful tools 
with which to search for dark matter. 


2.1 Neutrino Flux from Dark Matter Annihilation 

In general, the differential flux of monoenergetic neutrinos produced by dark matter annihi¬ 
lation can be expressed as 


d?®DM 

d9.dE 


47T 


-B v 



S(E U 


m x ), 


( 2 . 1 ) 


where dT/dV is the rate of dark matter annihilation per volume, and B u is the branching 
fraction for dark matter to annihilate to monoenergetic neutrinos, which we take to be flavor- 
independent. The differential antineutrino flux is the same. 

The Sun may be an important source of monoenergetic neutrinos from dark matter 
annihilations. If dark matter in the Sun is in equilibrium, then the annihilation rate is 1/2 

'in a liquid scintillation detector, the timing of illumination of the various photomultiplier tubes can be 
used to reconstruct the direction of the electron or positron track within the scintillator [10, 11]. 
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the rate at which dark matter is captured via scattering against solar nuclei (see, eg., [14]). 
We assume that dark matter-nucleon scattering is a spin-dependent contact interaction, as 
spin-independent scattering is tightly constrained by direct detection experiments. We may 
then express this capture rate [15] as T^ = Cq (rax)>< 0 g D x Kpx/pq)(v/270 krn/s) -1 ], where 
<7g D is the dark matter-proton spin-dependent scattering cross section, v is the dark matter 
velocity dispersion of a Maxwell-Boltzmann distribution, and where values of Cq D (rrix) may 
be found, for example, in [16]. We then find 
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where r® = 1 AU is the Earth-Sun distance, and the 5-function imposes the constraint that 
neutrinos will arrive from the direction of the Sun. If we take, as a benchmark, (gav) = 
1 pb x c and mx = 10 GeV, then the Sun is in equilibrium if cjg D > 3 X 10“ ' pb [17]. 
Comparing the result of Eq. 2.2 (assuming <jg D > 3 X 10“ ' pb) to the integrated annihilation 
rate near the Galactic Center as detailed in appendix A, it’s clear that the neutrino flux from 
the Sun is a few orders of magnitude larger than that from the Galactic Center * 2 . Henceforth, 
we will thus focus on searches for neutrinos arising from dark matter annihilation in the Sun 3 . 

The effects of oscillations, absorption, and regeneration as neutrinos propagate in the 
Sun, from the Sun to the Earth, and in the Earth, have been well-studied in recent years. 
Here, we consider only the energy range from 5 to 100 GeV, since evaporation from the 
Sun becomes significant for dark matter with mx 4 GeV [13] and neutrino absorption 
suppressed the line signal for neutrino energies E-, > 100 GeV [18]. For the range of energies 
we consider, although interactions in the Sun will affect neutrino flavor oscillations, the effect 
on the neutrino energy spectrum will be relatively small (see, for example, [19]). For a specific 
detector, a more detailed analysis could include the effect on the neutrino spectrum of matter 
effects in the Sun and Earth (the latter depend on the location of the detector), but for a 
general analysis this will not be necessary. 

Finally, for the range of WIMP masses considered here, we note that the assumption of 
dark matter capture-annihilation equilibrium is satisfied for total annihilation cross sections 
within a few orders of magnitude of the thermal benchmark, (gav) = lpb x c, for spin- 
dependent WIMP-proton scattering cross sections of a P SD ~ 10 -4 pb 4 . The total annihilation 
rate in the Sun is then set by g p sd as in Eq. 2.2, with a fraction B v of these annihilations 
resulting in monoenergetic neutrinos. Our results are therefore sensitive to the branching 
fraction to monoenergetic neutrino final states, but robust with respect to the overall dark 
matter annihilation cross section today, which may be significantly smaller than the typically- 
assumed thermal scale. 


2.2 Effective Area of the Detector 

Assuming that the number of protons and neutrons in the detector target are the same, the 
average neutrino-nucleon deep inelastic scattering cross sections for charged-current interac- 

“For mx 10 GeV, the neutrino flux from the Sun and from the Galactic Center both scale as mi“- 2 , so 
the conclusion holds for the range of dark matter masses considered here. 

4 Dark matter can only be captured in the Sun if it scatters against nuclei. However, this coupling need 

not imply that the dark matter annihilation branching fraction to hadronic final states is dominant. 

4 The equilibration time scales as (crg D < gav >) -1 ' 2 . 
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tions with E 3> GeV are [20] , 


<? V N ~ (6.66 x 10 3 pb) , 

Odn ~ (3.25 x 10- 3 pb) • (2-3) 


We may then express the effective area of the detector, A L JJ' = ov> ( ^ sL ) V, where pdet is 

ij v J\ y ” "iV y 

the density of the detector material, tun is the mass of a nucleon, and V is fiducial volume, 
for (anti)neutrino scattering as 


A e J f - = 2.03 x 10 _3 cm 2 
A% Sf ' = 9.88 x 10“ 4 cm 2 
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We will only consider events in which the 
where 


falls within an angle 0 con e of the Sun 


[ 21 ], 


e 


cone 


0.37 


10 GeV 
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(2.5) 


For electron neutrinos (antineutrinos) arriving from the sun, f u (/p) is the fraction of electrons 
(positrons) produced by charged-current interactions that fall within this cone. f U) p is roughly 
determined by the angular dependence of the matrix element for left-handed quarks to scatter 
against left-handed neutrinos and right-handed anti-neutrinos, respectively. f v p can then be 
expressed as 


fu 

fu 


l-p 2 r 1 1 
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( 2 . 6 ) 


where /3 is the boost from center-of-mass frame to lab frame. We find /„ ss 0.4 and /p ~ 0.8, 
results which are mildly dependent on the (anti)neutrino energy. In the following analysis, 
we use the full expression for given in eq. 2.6. 


2.3 Signal and Background Events 

The expected number of events in a particular detector arising from dark matter annihilation 
can then be expressed as 


Kun = 0.68 x T 




dE 


d®DM 

dE 



(2.7) 


where T is the detector runtime, max = mx[ 1 ± (e/2)], and e is the energy resolution of 
the detector. Again, B u is the branching fraction for dark matter annihilation to all neutrino 
flavors; since we assume that annihilation is flavor-independent, ~ 1/3 of the neutrinos 
reaching the detector will be electron neutrinos even after including vacuum and matter 
oscillation effects [22], The factor of 0.68 accounts for the fact that 68% of events fall within 
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a bin of width e. Similarly, the number of events due to the atmospheric neutrino background 
can be expressed as 


KtZn = T 


dE 


d 2 $ 


atm 


d9dE 


x A e J/' x [(27 t)( 1 - cos 0 cone )\ 


( 2 . 8 ) 


where d d ^ d g is the atmospheric electron (anti)neutrino flux. We use the angle averaged 
atmospheric electron (anti)neutrino fluxes given in [23], which, in the relevant energy range, 
are well-approximated by the power law fits 


d 2 Ktm 

d9.dE 

d 2 €tm 

d9dE 


(4.17 x 10~ 2 cm _2 s” 1 sr _1 GeV” 1 ) x (°- 80 + 
(2.42 x 10 _2 cm~ 2 s“ 1 sr _1 GeV _1 ) x ^0.53 + 


-3.490 


-3.417 


(2.9) 


Given the number of signal and background events expected at a neutrino detector, one 
can then determine sensitivity of the detector to a dark matter model for any choices of mx 
and <Tg D . 


3 Bounds from KamLAND 

In this section, we apply the analysis framework above to the specific case of the KamLAND 
detector and the data it has accumulated over more than 10 years. KamLAND [24] is a 
liquid scintillation neutrino detector with p de t ~ 0.8 g/cm 3 and an approximately spherical 
inner detector with volume ~ 1000 m 3 . We define the fiducial volume of the detector such 
that if a charged-current interaction produces an within the fiducial volume, at least 10 
radiation lengths (~ 4.3 m) will be contained within the inner detector. For KamLAND, the 
fiducial volume is ~ 500 m 3 [17]. If an electron (anti)neutrino with energy greater than 1 GeV 
participates in a charged-current interaction within the fiducial volume, it is estimated that 
the energy of the original neutrino can be reconstructed to within a few percent accuracy, 
and that approximately 3600 live-days of data are available for such an analysis [25]. As 
a benchmark, we will take the energy resolution for this analysis to be e = 5%. It is also 
estimated that the charged lepton flavor can be determined with very high efficiency, and 
that the direction of the electromagnetic shower can be determined with a resolution much 
better than 0 CO ne for the energy range of interest [11]. 

We note that for the exposure of KamLAND, one expects far less than one background 
event within each energy bin, whose size is determined by the energy resolution. In Figure 1, 
we show the expected number of background events as a function of the central bin energy 
for a KamLAND exposure of 3600 days and energy resolution of 1% (solid grey), 5% (solid 
black), and 10% (dotted grey). We also show the projected number of background events as 
a function of energy for a detector with an exposure 10 times larger than that of KamLAND 
with an energy resolution of 3% (dashed black). The inset shows the low-energy regime 
in more detail. It is clear that, for the exposure of KamLAND, the atmospheric neutrino 
background is negligible. 

A 90% confidence level upper limit on the spin-dependent dark matter-proton elastic 
scattering cross section as a function of dark matter mass is obtained (see, for example, 
[26]), assuming zero observed events in each relevant energy bin during 3600 KamLAND 
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Figure 1. The expected number of background events as a function of the central bin energy for a 
KamLAND exposure of 3600 days with energy resolution of 1% (solid grey), 5% (solid black), and 
10% (dotted grey). We also show the projected number of background events as a function of energy 
for a detector with an exposure 10 times larger than that of KamLAND with an energy resolution 
of 3% (dashed black). We note that the dotted grey curve also represents the expected number of 
background events for a detector with an exposure 10 times larger if the energy resolution is 1%. 


live-days. The expected upper limit on the spin-dependent dark matter-proton scattering 
cross section from KamLAND is shown in Figure 2 as the thick black contour for B u = 1, 
px = Pq and v = 270 km/s 5 , as well as the expected sensitivity of a liquid scintillation 
detector with a factor 10 larger exposure than KamLAND and with a 3% energy resolution 
(thick black dashed). For comparison, the dashed, solid, and dotted magenta curves are the 
upper limits from Baikal NT200 detector [27], assuming dark matter annihilation to F e z/ e , 
and z/ T z/ T , respectively, the solid cyan curve is the upper limit from Super-Kamiokande, 
assuming annihilation to t + t~ [28], and the dashed blue contour is the upper limit from 
the PICO-2L Bubble Chamber dark matter experiment [29]. For clarity, we display only 
the leading published results (or, in the case of PICO-2L, publicly released but not yet 
published), omitting competitive but subleading bounds from direct and indirect searches 
such as PICASSO [30], SIMPLE [31], COUPP [32], the Baksan Underground Scintillator 
Telescope [33], and the IceCube Neutrino Telescope [34]. 

The estimated sensitivity is relatively robust, even if the actual energy resolution is 
different from our benchmark estimate by a factor of a few. This is clear from Fig. 1 for 
the current KamLAND exposure of ~ 3600 days. Even for an exposure 10 times larger, the 
atmospheric background only approaches 0(1) event in an energy bin for the very lowest 
energies considered. For e = 3%, the expected number of background events in an energy 
bin is never larger than 0.39, for E u = 5 GeV, and falls to less than 0.09 for E u > 10 GeV. If 
a single event is observed in an energy bin, the actual sensitivity will drop by ~ 40 — 50%. 

5 As the number of signal events is proportional to ag D -B„(pA'/p©)(h/270 km/s) _1 , the sensitivity for any 
choice of B v , px and v can be determined by a simple rescaling. 
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Figure 2. 90% CL upper limit on the spin-dependent dark matter-proton elastic scattering cross 
section as a function of dark matter mass assuming zero events at KamLAND in 3600 days of data 
(thick black), as well as the expected sensitivity of a liquid scintillation detector with a factor 10 larger 
exposure than KamLAND and with a 3% energy resolution (thick black dashed). We assume B u = 1, 
Px = Po and v = 270 km/s. For comparison, the dashed, solid, and dotted magenta curves are 
the upper limits from Baikal NT200 detector [27], assuming dark matter annihilation to 
and v T v T , respectively, the solid cyan curve is the upper limit from Super-Kamiokande, assuming 
annihilation to t + t~ [28], and the dashed blue contour is the upper limit from the PICO-2L Bubble 
Chamber dark matter experiment [29]. 


As is clear from Figure 2, at this moment, the limit obtainable from KamLAND data 
would be stronger than limits from traditional direct dark matter searches (i.e. for dark matter 
collisions with nuclei in terrestrial detectors) for WIMP masses <15 GeV, assuming dark 
matter annihilates to vv final states with B u = 1. The sensitivities of neutrino telescopes such 
as Super-Kamiokande, Baksan, and IceCube are specific to the dark matter annihilation final 
state, though harder neutrino spectra generally yield stronger constraints on the scattering 
cross section. For those detectors, actual limits for monoenergetic neutrino final states are 
likely stronger than the limits shown here. We encourage collaborations to perform these 
searches. 

Finally, we note that KamLAND’s sensitivity is signal-limited, due to its relatively 
small exposure compared to that of much larger water Cherenkov detectors such as Super- 
Kamiokande. As such, the utility of liquid scintillation detectors for searches for dark matter 
annihilation to monoenergetic neutrinos can only be fully realized for detectors with a much 
larger exposure. For example, for an exposure ten times larger than our benchmark Kam¬ 
LAND exposure, cross sections as small as a p SD ~ 10 -4 (6 X 10~ 4 ) pb could be probed for 
mx = 10 (100) GeV. Such an experiment would be an very powerful probe of the interactions 
of dark matter with nuclei in the Sun. This is especially true at low WIMP masses, where 
liquid scintillation detectors have robust sensitivity to the WIMP-proton scattering cross 
section down to mx ~ 4 GeV, below which evaporation from the Sun becomes significant. 
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4 Conclusions 


In this work, we have considered the sensitivity of liquid scintillation neutrino detectors to 
dark matter annihilation in the Sun in the case that the dark matter annihilates to rnonoen- 
ergetic neutrinos. Because liquid scintillation detectors can reconstruct with good accuracy 
all of the energy of an electron (anti)neutrino interacting through a charged-current interac¬ 
tion, this is essentially a search for an excess in a single energy bin. The resulting reduction 
in background greatly enhances the sensitivity of liquid scintillation neutrino detectors. As 
an example, we have considered the sensitivity that KamLAND may be capable of, using 
existing data. 

As may be expected, the sensitivity of liquid scintillation detectors found here exceeds 
those of other estimates [35], which focussed on annihilation to final states that included a 
continuum spectrum of neutrinos, in which case a signal would emerge over a broad range 
of neutrino energies rather than within a single energy bin. With its current exposure, 
KamLAND would expect far less than one background event in any single energy bin for 
energies > 5 GeV. In fact, for any liquid scintillation detector, the expected sensitivity will 
increase linearly with exposure in this low-background regime (this low-background regime 
should extend up to an exposure of roughly 10 times the current exposure of KamLAND for 
all WIMP masses considered, and to considerably larger exposure for rrix > 100 GeV). If the 
exposure becomes large enough that the expected number of background events in an energy 
bin approaches one, the assumption of zero events in any bin is, of course, no longer valid 
for estimating the sensitivity. 

It is worth noting that liquid argon-based neutrino detectors can also reconstruct the 
energy and direction of an electron shower created by a charged-current interaction. A 
similar analysis could thus be performed for such detectors. Moreover, future liquid argon- 
based detectors may have a much larger exposure than KamLAND. For example, the LBNE 
far detector is foreseen to be a 34 kiloton liquid argon time projection chamber [12]. If the 
energy resolution is as good as a few percent, then such detectors could obtain competitive 
sensitivity, placing them among the leading experimental searches for particle dark matter. 
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A Monoenergetic Neutrinos from the Galactic Center 

Indirect signals of dark matter annihilation in most astrophysical bodies do not rely on dark 
matter scattering with nuclei in order to become gravitationally bound, as is the case for 
the annihilation signal from the Sun. Instead, the annihilation rate is proportional to the 
annihilation cross section, rather than the dark matter-nucleon scattering cross section. For 
dark matter annihilation in the Galactic Center region, the integrated annihilation rate may 
be expressed as 
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where (ctav) is the total dark matter annihilation cross section times velocity, rg = 8.33 kpc 
is the distance from the earth to the Galactic Center [36], and pg = 0.3 GeV/cm 3 is the 
putative density of dark matter in the solar system [37]. The quantities r & and pg serve to 
normalize the J-factor, which encodes the relevant information regarding the dark matter 
density along the line of sight to the target. Typical J-factors for the Galactic Center are 
0(1 — 10) [18]. We have assumed that the dark matter particle is its own anti-particle; if 
the particle and anti-particle are distinct, this annihilation rate is reduced by a factor of 
2. The differential flux of monoenergetic neutrinos produced by dark matter annihilation 
in the Galactic Center (or other region, given the appropriate J-factor) may be calculated 
according to equation 2.1 with the integrated annihilation rate as in A.l. 
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